What is wrong with the Lax - Richtmyer Fundamental Theorem of Linear Numerical Analysis ?
نویسنده
چکیده
We show that the celebrated 1956 Lax-Richtmyer linear theorem in Numerical Analysis-often called the Fundamental Theorem of Numerical Analysis-is in fact wrong. Here " wrong " does not mean that its statement is false mathematically, but that it has a limited practical relevance as it misrepresents what actually goes on in the numerical analysis of partial differential equations. Namely, the assumptions used in that theorem are excessive to the extent of being unrealistic from practical point of view. The two facts which the mentioned theorem gets wrong from practical point of view are :-the relationship between the convergence and stability of numerical methods for linear PDEs,-the effect of the propagation of round-off errors in such numerical methods. The mentioned theorem leads to a result for PDEs which is unrealis-tically better than the well known best possible similar result in the numerical analysis of ODEs. Strangely enough, this fact seems not to be known well enough in the literature. Once one becomes aware of the above, new avenues of both practical and theoretical interest can open up in the numerical analysis of PDEs.
منابع مشابه
Ju l 2 00 5 What is wrong with the Lax - Richtmyer Fundamental Theorem of Linear Numerical Analysis ?
We show that the celebrated 1956 Lax-Richtmyer linear theorem in Numerical Analysis-often called the Fundamental Theorem of Numerical Analysis-is in fact wrong. Here " wrong " does not mean that its statement is false mathematically, but that it has a limited practical relevance as it misrepresents what actually goes on in the numerical analysis of partial differential equations. Namely, the as...
متن کاملApplication of Convolution of Daubechies Wavelet in Solving 3D Microscale DPL Problem
In this work, the triple convolution of Daubechies wavelet is used to solve the three dimensional (3D) microscale Dual Phase Lag (DPL) problem. Also, numerical solution of 3D time-dependent initial-boundary value problems of a microscopic heat equation is presented. To generate a 3D wavelet we used the triple convolution of a one dimensional wavelet. Using convolution we get a scaling function ...
متن کاملConsistency & Numerical Smoothing ⇒ Error Estimation — An Alternative of the Lax - Richtmyer Theorem
We all know Lax-Richtmyer Theorem: Assuming consistency, convergence is equivalent to numerical stability. However, it is in practice very difficult to verify the numerical stability of a scheme while solving an evolution equation, especially if the equation is nonlinear and/or the scheme is complex. Consequently, a large gap exists between error analysis theory and numerical computation practi...
متن کاملA Partial Nonlinear Extension of Lax-Richtmyer Approximation Theory
A partial nonlinear extension of Lax-Richtmyer approximation theory by Aradhana Kumari Adviser: Professor Dennis Sullivan Lax and Richtmyer developed a theory of algorithms for linear initial value problems that guarantees, under certain circumstances, the convergence to numerical solution of initial value problem. The assumptions are first that the difference equations (algorithms) approximate...
متن کاملQuantitative Theory of Richtmyer-meshkov Instability in Three Dimensions
A material interface between two fluids of different density accelerated by a shock wave is unstable. This instability is known as Richtmyer-Meshkov (RM) instability. Previous theoretical and numerical studies primarily focused on fluids in two dimensions. In this paper, we present the studies of RichtmyerMeshkov instability in three dimensions in rectangular coordinates. There are three main r...
متن کامل